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Abstract 

5— i ■ In this Letter we discuss a natural general relativistic mechanism 

^Q. that causes inhomogeneities and hence generates matter perturbations 

in the early universe. We concentrate on spikes, both incomplete 
spikes and recurring spikes, that naturally occur in the initial oscilla- 
{SJ ■ tory regime of general cosmological models. In particular, we explicitly 

show that spikes occurring in a class of G2 models lead to inhomo- 
geneities that, due to gravitational instability, leave small residual 
imprints on matter in the form of matter perturbations. The resid- 
^5 \ ual matter overdensities from recurring spikes are not local but form 

CN on surfaces. We discuss the potential physical consequences of the 

residual matter imprints and their possible effect on the subsequent 
formation of large scale structure. 



> 



(N 



X 



Introduction. It is a general feature of solutions of partial differential 
equations (PDE) that spikes 1 occur [1]. Therefore, spikes are expected to 
occur in generic solutions of the Einstein field eqns (EFE) of general relativity 
(GR), which are a complicated set of PDE [2]. Indeed, when a solution of 
the EFE is stable at certain points but otherwise unstable, spikes will arise 
near the stable points. The set of such points can be a surface, a curve or a 
point in the three-dimensional space. 

Belinskii, Khalatnikov and Lifshitz (BKL) [3] have conjectured that within 
GR, the approach to the generic (past) spacelike singularity is vacuum dom- 
inated, local, and oscillatory (i.e., Mixmaster) . Studies of G2 2 and more 
general cosmological models have produced numerical evidence that the BKL 
conjecture generally holds except possibly at isolated points (surfaces in the 



1 Spikes are narrow inhomogeneous structures (see, e.g., Figure 2). 
2 The so-called G2 spacetimes are those which admit two commuting spacelike Killing 
vector fields, and hence only have one degree of inhomogeneity. 
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three-dimensional space) where spiky structures ('spikes') form [4]. These 
spikes become ever narrower as the singularity is approached. The presence 
of such spikes violates the local part of the BKL conjecture. 

The study of spikes is severely limited due to the enormous numerical re- 
sources needed to resolve the narrowing spikes in simulations. In [5], further 
improved numerical evidence was presented that spikes in the Mixmaster 
regime of G<i cosmologies are transient and recurring, supporting the con- 
jecture that the generalized Mixmaster behavior is asymptotically non-local 
where spikes occur. It is believed that this recurring violation of BKL local- 
ity holds in more general spacetimes; however, it remains to study this more 
comprehensively. 

In this Letter, we wish to study the residual imprints of the spikes on 
matter inhomogeneities. As the spike inhomogeneities form, matter under- 
goes gravitational instability and begins to collapse to form overdensities. 
In this way, density inhomogeneities naturally form in the early universe in 
generic GR models. This is a natural GR effect, and the question arises as 
to whether these density inhomogeneities are physical and can be observed 
and, indeed, whether such matter inhomogeneities could act as seeds for the 
subsequent formation of large scale structure. 3 

Therefore, we are interested in possible (purely classical GR) effects that 
could cause inhomogeneities and hence generate matter perturbations that 
might then impinge on structure formation. We are particularly interested 
in recurring and distributed spikes formed in the oscillatory regime (or re- 
curring spikes for short), and their imprint on matter. We show that matter 
inhomogeneities can occur in simple so-called orthogonally transitive (OT) 
G 2 models [4] with a tilted radiation fluid. 4 Some intuition can be gained 
by looking at the case of small (normalized) density parameter Q with small 
tilt. First, we study a test fluid in an exact vacuum background [6] (small 
tilt and small Q with no backreaction on the geometry); although a complete 
spike transition leaves no imprint in this partial spike transition does 

leave an imprint. We then study the small tilt approximation with small Q 
(when matter is not a test field), and show heuristically that residual matter 
inhomogeneities occur. The argument doesn't apply to the typical case of 
non-negligible Q and/or tilt; however, we argue it does occur in more general 
models (partly by studying the decay rates of the inhomogeneous spikes). 

3 We shall use the general terminology matter inhomogeneities or density perturbations 
or overdensities to describe the residual local inhomogenieties imprinted in the matter, 
but the actual mathematical behaviour is rather more complicated and perhaps can be 
better descibed as a delay in evolution at the spike (also note that fi and the physical fluid 
density p may have different qualitative behaviour, due to the evolution of the Hubble 
normalization). 

4 A tilted fluid is a fluid that has a non-zero velocity (or tilt) relative to the chosen 
reference frame. 
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We also discuss some simple exact spike solutions (previously studied in the 
literature), which have residual matter inhomogeneities, further lending sup- 
port to the heuristic analysis [2]. But ultimately we shall demonstrate the 
existence of matter perturbations by numerical simulations. 

Both the incomplete spikes and the recurring spikes are potentially of 
physical importance (but perhaps for different reasons). We are particularly 
interested in the (G2) recurring spikes, and we show explicitly that there 
exist spikes leading to inhomogeneities and a small residual in the form of 
matter perturbations. Moreover, these density perturbations occur naturally 
within generic cosmology models within GR and they are not local (or point- 
like) but form on surfaces and give rise to a distribution of perturbations. 
These classical GR generated inhomogeneities consequently may affect large 
scale structure. Larger effects (e.g., occuring from incomplete spikes), must 
of course be consistent with observations. 

In the final section we discuss the potential physical applications of the 
results, and speculate on whether these recurring spikes might be an alter- 
native to the inflationary mechanism for generating matter perturbations. 

Spike Analysis. The full evolution (EFE) equations are given in Ap- 
pendix D of [2] (with A = 0). The variables there are \/3\ -normalized (or 
Hubble-normalized) . The (5- normalized lapse is chosen to have the value 
— 1/2 here (which is appropriate for numerical simulation). The evolution 
equations are: 5 

d T \n\$\ = ^[l + S 2 _ + S 2 x +iV^+iV! + ( 7 -l)fi] (1) 
d T In Q = -^vE-^dx In Q + -'yE^dxV 

-?(2 - 7 )[1 + S 2 + S 2 + Nl + Nl - fl], (2) 

where S + and q (and all terms) are given in [2]. There are also evolution 
equations for Ei l and v (and E_, £ x , iV_, N x ). In this Letter we need only 
focus on Q and v. We consider the linearization of each variable around a 
background metric: 

n = fio + e^i + 0(e 2 ), v = v + ev 1 + 0(e 2 ) (3) 

(etc.). In the small Q and small v approximation we assume that Q and v 
have vanishing zeroth order terms. The linearized evolution equations are 
easily obtained from the full EFE. 

5 f3 is the area expansion rate of the (y, z) plane; 7 is the equation of state parame- 
ter, with 7=5 describing the radiation fluid. S + , £_ and E x are components of the 
/3-normalized rate of shear; N x and N_ are components of the /3-normalized spatial cur- 
vature, is the /3-normalized fluid density; v is the relative fluid velocity (tilt) in the 
x-direction. /3-normalization is analogous to the standard Hubble-normalization, and is 
related through H = — E+). 
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The expressions for the exact vacuum spike solution [6], which are used 
as the zeroth order (background) solution in the linearization are: 

(£_, N x , £ x , AL) = ^-cS_ Taub - -J=, sAL Taub , cAL Taub , -sS_ Taub ^ . (4) 



r - 1 2/ 



/ 2 + l' / 2 + l' 



/ = we r sech(wr)x. (5) 



w 1 w 

S-Taub = -^=tanh(wr) - -^=, AL Tau b = -^=sech(wr) (6) 

More importantly, we note that 

(5 = -isech( W r)e^-^ A2 (/ 2 + l)~i (7) 

We shall provide analytical evidence for the existence of residual mat- 
ter inhomogeneities due to spikes. We study the case of negligible Hubble- 
normalized density f2 and negligible tilt v. To see the effect of a spike on 
matter, we first consider a test perfect fluid. 

Spike imprint: test fluid and incomplete spikes: For a test fluid (with 
negligible Q and v), from the EFE and the conservation equation we obtain: 

d T ln\(3\ = ^q + l) (8) 

d T lnfi = -(g + l) + ^ 7 (l-£ + ). (9) 
Remarkably, (p and) f2 can be solved exactly as follows: 

'='° w G4) 7 - a ^ {x) (mT 2 '"- (io) 

Note that d T In \/3\ is positive. So for 7 satisfying < 7 < 2, we see that 
(3 and p blow up at the singularity, but Q = p/(3/3 2 ) tends to zero ["matter 
does not matter"]. As spacetime expands, p and /3 decrease and Q increases. 

For the spike solution, j3 is given by (7), where 

(f + = [(we T sech(wr)x) 2 + l]"*, (11) 

depends on x and so it is responsible for the spatial inhomogeneity. It can be 
seen from Fig. 1 that this factor is asymptotically homogeneous (and equals 
1) as r tends to ±00, for \w\ > 1. It takes the value 1 at x = for all time, 
but is smaller for x^O. This means that f3 is inhomogeneous during a spike 
transition, but a complete spike transition restores j3 to homogeneous. So 
the cumulative effect of a complete spike transition on spatial inhomogeneity 
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Figure 1: The factor (f 2 + with w = 3. x and r are the dimensionless 
space and time variables, respectively, r tends to +oo at the singularity. The 
spatial axis used here (and in Figures 2 and 3) is X = e T x. 



in f3 is zero. Similarly, a complete spike transition has zero cumulative effect 
on the inhomogeneity in p and fl (when regarded as a test fluid). 

Discussion: It is interesting that for a partial spike transition (i.e., during 
a spike transition), (3 develops a hump at the spike surfaces x — 0. Corre- 
spondingly, p develops a hump at x — 0, and fl develops a dip at x — 0. 
This can be interpreted as a delay or lag in the evolution at the spike point 
compared with points far away. The delay is temporary, and by the end of 
the spike transition, there is no delay. 

As the universe expands from the big bang, presumably it starts with 
a Mixmaster oscillatory regime with spike transitions occuring on surfaces 
here and there, and with Bianchi type II transitions occuring elsewhere. 
Eventually, the oscillatory regime ends when fl is no longer negligible. Some 
of the spike transitions are in the middle of transitioning when the oscillatory 
regime ends. This leaves an inhomogeneous imprint on matter as well as 
curvature, by way of a delay in evolution at surfaces of spike points. and 
p are larger, and fl is smaller at surfaces of spike points than at other points. 

Since fl is no longer negligible at the formation of an incomplete spike, 
eqn. (10) may not be a good approximation. But the conclusion that \(3\ 
and p are larger, and fl is smaller at surfaces of spike points, is confirmed 
by numerical simulations. We note that the isolated incomplete spikes are 
rare, since the solutions spend much less time undergoing transitions than 
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Figure 2: \(3\ at r = with w = 3, A 2 = 0. 



evolving in the Kasner epoch. 

Heuristic analysis: small Q. We next study the linearized equations with 
a spike background, assuming a small Q. The zeroth order terms in the 
linearized equations are satisfied identically by the exact spike solution. As- 
suming Q Q = 0, then the leading order term in the linearized EFE gives 
in terms of (3q (eqn. (7)): 



where C is a constant. This means the initial spatial profile of Qi is not 
freely specifiable (for non-dust), but is determined by the spike background. 

We can then solve for (Ei l )i. The other unknowns, namely (3\, (£-)i, 
(S x )i, (AL)i and (iV x )i, are coupled, although v\ decouples from these un- 
knowns. However, even in this case, where we assume that both Q and v are 
small, we must ultimately resort to numerical analysis. 

Heuristic analysis: non-trivial Q. Let us now consider the large Q case 
with Q 0. From these eqns (for vq small) we obtain 




(12) 




d T \nQ = 8 T In |/3|- 2 " 7 + ^ 7 (2 - 7 )0. 
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Figure 3: Integral of B; see Fig. 1. 



Writing Q = fi (l + efti), where Q = Q (x)f3 Q 
solution as given by eqns. (10) and (7), we obtain 



(2-7) 



is the zeroth order 



and hence 



i 7(2- 7 )fi (l + efi 1 ) 2 , 



37(2-7) 



~i -1 



Q cLt 



where Q = Ci (x)F{t)B(t,x), and B(t,x) = (1 + f 2 f-^/ 2 (and F(t) is 
defined by eqn. (7) and / is defined by (11)), plus possible additional con- 
tributions from 0x. The important point is that there are consequently con- 
tributions to Qi from terms like J B(r,x)dr. Now, the spike occurs for 
— r < r < r , but as r — > 00 the spike disappears and B{r,x) becomes 
homogeneous (although transient inhomogeneities occur for — r < r < r , 
see Fig. 1) However, terms like j B(r,x)dr retain a residual inhomogeneity. 
Hence terms like J B(t, x)dr (and integrals thereof) contribute an inhomo- 
geneous imprint to Q, and hence the matter density. This is illustrated in 
Fig. 3 (for the same B as in Fig. 1). 

Therefore, there will be a residual inhomogeneous imprint on the density 
due to a spike. This has been illustrated here in this simple case, but it is 
most likely to be true in all generality. We still need to analyse the case of 
large Q and large v. However, this can only be done numerically. 

Decay rate: From above Qq = Q (x)F(t)B(t,x), where 



B(t,x) = [1 + (we T sech(wr)x 



,21 1-7/2 



(7 < 2). The residual inhomogeneity comes from B(t,x). Of course, as 
t — > 00, B(t,x) — > 1, and the inhomogeneity from the spike decays. For 
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large r (or equivalently for r — > — oo), we have that 

B(t, x) ~ 1 + 4w 2 (l - 2) e 2(1 - w)r a: 2 

where w > 1 (assuming w positive). Hence the inhomogeneity (the imprint 
of the spike) decays at the rate e 2 ^ 1 ~ w ^ T . 

Now, from linear perturbation theory, density inhomogeneties are ex- 
pected to grow at the rate t c ~ e CT , for some constant c, due to gravitational 
instability (where c depends on 7; for example, and simply as a reference, in 
a flat Friedmann-Lemaitre (FL) model, c = (7 — |) on the largest scales). 
Thus, it might be expected that residual inhomogeneous imprints on the 
density would occur whenever c + 2(1 — w) > (i.e., for 1 < w < 1 + c/2). 

Exact solutions with matter; LRS G2 dust and LTB models: Finally, we 
briefly discuss some simple exact spike solutions that have been studied pre- 
viously, which constitute examples of residuals in special cases. 

Spikes arise when a solution straddles the stable manifold (or separatrix) 
of an unstable equilibrium point (a source or a saddle point) [2, Chapter 6]. 
If the unstable manifold of that equilibrium point is one-dimensional, then 
spikes occur on a surface (in three-dimensional space); if two-dimensional, 
then on a curve; if three-dimensional, then on a point. 

In [2], an explicit solution in the class of locally-rotationally symmetric 
(LRS) G 2 cosmological models with dust and a cosmological constant (a 
special case of the Szekeres solution) was given, illustrating explicitly how 
the straddling of the stable manifold of a saddle point leads to the formation 
of a spike. Solutions on each side of the stable manifold approach different 
sinks, leading to a discontinuous limit in some of the variables. The shape 
of this discontinuous limit may look like a step function or a spike. The 
explicit form for the inhomogeneous Q(x) (with A = 0), which develops a 
step/jump inhomogeneity, was given in [2]. Similarly, it is possible to find 
a spiky solution in Lemaitre-Tolman-Bondi (LTB) models [2], which may be 
more relevant to structure formations than LRS G2 models (we shall persue 
this further in future work). 

Numerical evidence of residuals. Let us do a simple numerical simu- 
lation to show that there is indeed an inhomogeneous imprint left by a spike. 
The numerical code is essentially the one used in [5], with a radiation fluid 
added. The limitation of the code is that it cannot handle shock waves and 
the step-like structure that forms in the tilt. The tilt is unstable for S + < 0, 
so this limits the simulations to spikes with roughly \w\ < 2. The current 
code with zooming also means that each simulation can only see the space- 
time at horizon scale. To see the imprint, which becomes super-horizon into 
the past, several simulations are needed to produce the figure. It is inef- 
ficient, and we hope to find an efficient way to simulate the spacetime at 
super-horizon scale in the future. 
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Figure 4: Ratio of f2 along different worldlines. See text for details. 

Fig. 4 shows the time evolution of 7 separate simulations centred on 7 
worldines 



(labeled worldline number 6 to in the plot). The ratio of Q is taken to be 



Hence the plot along woldline number equals 1. From the plot, towards the 
singularity we see that Q near the spike (at x = 0) grows 6-fold. Reversing 
the time direction, we see that the spike thus leaves an imprint on Q with a 
6-fold underdensity. The above plot is produced with Q = 1CT 10 and w — 1.5. 
Using w = 1.4 yields a 5.8-fold change. Using Q = 1CU 5 does not affect the 
result. From these preliminary results, we conclude that the spike leaves an 
underdense imprint on Q regardless of the size of Q. The amplitude of the 
underdensity depends on w. 

Let us make two comments about the imprint obtained. First, the wave- 
length is at the width of the spike (roughly the horizon scale) when it is 
created. But towards the future it will become sub-horizon, as the horizon 
expands. Because of this, it is difficult to see the imprint numerically when 
using a single simulation with zooming. Second, the imprint is large (6-folds) 
because Q is close to zero. Towards the future this ratio will become smaller 
as Q becomes of order 1 in the radiation dominated close-to-flat-FL era. 



x = 0, 1(T 7 , 1(T 6 , 1(T 5 , 1(T 4 , 1(T 3 and 10 -2 



(14) 




(15) 
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We conclude that the numerics show the occurance of residuals (i.e., def- 
inite spatial dependence is illustrated). We are not too concerned at this 
time about the shape/characteristics of the residuals. However, we can ask 
whether the numerics suggest void formation (see Fig. 4). Q develops a void 
at a spike location; they are voids when they form. But it is hard to say 
what these imprints will lead to, since they may change to overdensities via 
subsequent later dynamics. 

Discussion. We are interested in the possible existence of a GR mecha- 
nism for generating matter perturbations. We have concentrated on spikes, 
both incomplete spikes and recurring spikes, and shown that there are ef- 
fects, entirely within classical GR, that could cause inhomogeneities and 
hence matter perturbations. 

In particular, we have shown that there will be residual matter perturba- 
tions from spikes, based on a heuristic qualitative analysis of a single spike in 
the OT G 2 model and exact LRS G 2 and LTB spike solutions with matter [2], 
and most importantly from numerical simulations. There are residuals from 
an incomplete spike, that might in principle be large and thus affect structure 
formation (and any such effects might lead to observational constraints). 

In addition, we have explicitly shown that there exist G 2 recurring spikes 
that lead to inhomogeneities and a residual in the form of matter pertur- 
bations, that these occur naturally within generic cosmology models within 
GR, and that they are not local but form on surfaces and give rise to a distri- 
bution of perturbations. In the G 2 models the inhomogeneities can occur on 
a surface, and in general spacetimes the inhomogeneities can occur along a 
line, leading to matter inhomogeneities forming on walls or surfaces. Indeed, 
there are tantalising hints that (from dynamical and numerical analyses) that 
filamentary structures and voids would occur naturally in this scenario. 

Also there are shock waves that could generically form (generated by 
inhomogeneities in the fluid pressure and density, even within inhomogeneous 
Newtonian cosmology), not necessarily associated with spikes, and leave a 
residual imprint in the matter perturbations. Incomplete spikes, or shock 
waves not associated with spikes, might leave isolated (and perhaps large) 
perturbations. However, shocks might form associated with spikes (indeed, 
spikes themselves will likely cause shock waves which will, in turn, give rise 
to additional inhomogeneities and hence matter perturbations), and there 
might be interesting behaviour due to interaction of a spike and a shock 
wave forming close to a spike, leaving additional residual perturbations [2]. 

Inflationary cosmology provides a causal mechanism which generates the 
primordial perturbations which were later responsible for the formation of 
stars, galaxies, clusters, and all large scale structures of our Universe un- 
der the influence of gravitational collapse. The density perturbations pro- 
duced during inflation are due to quantum fluctuations in the matter and 
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gravitational fields [7]. The evolution outside the Hubble radius then pro- 
duces a large amplification of the perturbations. Primeval fluctuations are 
then thought to be present at the end of the inflationary epoch. Provided 
that inflation lasts sufficiently long, generically an (almost) scale-invariant 
Harrison-Zel'dovich spectrum of density fluctuations are generated, which 
then evolve to the (tiny, 10~ 5 ) adiabatic, Gaussian, and scale invariant den- 
sity fluctuations in the power spectrum of the CMB [7] (which then allows 
for a detailed comparison with current observations). 

If inflation does occur after the mixmaster regime is over, all of the clas- 
sical inhomogeneities will be redshifted away, but any 'intermediate' scale 
effects might not necessarily be redshifted away, and might still be of impor- 
tance in cosmology. In the standard scenario it is assumed that there is no 
scale between the quantum scale and the classical scale, that would be am- 
plified during inflation and become important in cosmology after inflation. 
We could ask whether there is such an intermediate scale inherent in GR. 

Perhaps the scale associated with distributed recurring spikes might be 
such an intermediate scale. Indeed, since the 'horizon' goes to zero as we 
approach the initial spacelike mixmaster 'singularity' in the classical regime, 
the scale (wavelength) of the 'first' recurring spikes can be arbitrarily small 
(but larger than the Planck scale assuming the mixmaster oscillations occur 
in the classical regime after leaving the Planck regime) . If there are classical 
scale inhomogeneities produced (either by late time, isolated spikes or shocks, 
or amplified recurring spikes) they must be consistent with current CMB 
observations. 

Indeed, in the case of recurring spikes within GR, the perturbations can 
be extremely small (and, as noted earlier, occur in non-isolated distributions 
everywhere). Therefore, we speculate whether recurring spikes could generate 
primordial matter perturbations that are very small (10~ 5 ) at recombination 
and subsequently seed the large-scale structure of the actual universe, and 
consequently act as an alternative to the usual inflationary mechanism. In 
spite of the remarkable success of the inflationary Universe paradigm, there 
are several serious conceptual problems for current models [8]. Therefore, a 
classical GR mechanism for generating small inhomogeneities would be an 
interesting alternative, and might perhaps bring GR back to the centre of 
cosmology [9]. This is, of course, highly speculative, but it might motivate 
the further study of spikes in GR. 

Therefore, in future work we shall study what happens in more generality, 
hopefully obtaining some information regarding the properties of residual 
matter inhomogeneities (e.g., their characteristic scales) and the statistical 
properties of matter perturbations (i.e., the distribution of perturbations 
in general inhomogeneous models). In G2 models perturbations occur on 
planar surfaces, while in general inhomogeneous (G ) models perturbations 
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may occur on other types of surfaces. 

Unfortunately, although this is very interesting, this is also extremely 
difficult to investigate, since it involves the numerical integration of general 
GR cosmological models. The numerical challenges include a need to develop 
a numerical code capable of simulating spacetimes at super-horizon scales 
while at the same time able to resolve spikes at sub-horizon scales, and to 
develop a similar numerical code capable of handling shock waves. 

Acknowledgment: This work was supported, in part, by NSERC of Canada. 
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